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A  CALCULATION  OF  DESCENT  SEGMENT  OF  INTERNAL  BALLISTICS  IN 
SOLID  ROCKET  MOTORS 

Fang  Dingyou 

Abstract 

The  precise  prediction  of  the  after-working  impulse  of 
solid  rocket  motors  is  very  important  for  the  prediction  of 
missile  and  motor  performance.  This  paper  uses  the  character¬ 
istic  method  to  consider  the  changes  of  gaseous  parameters 
along  the  combustion  chamber  length  and  to  calculate  the 
descent  segment  of  internal  ballistics  in  solid  rocket  motors. 
The  gaseous  parameters  in  the  descent  segment  along  the  com¬ 
bustion  chamber  length  which  change  with  time  are  given  and 
the  calculation  results  of  the  characteristic  method  and 
isentropic  model  are  compared. 

I.  Preface 

The  precise  prediction  of  the  after-working  impulses  of 
solid  rocket  motors  is  very  important  for  the  prediction  of 
missile  and  motor  performance.  In  the  past,  calculations  of 
the  descent  segments  of  internal  ballistics  in  solid  rocket 
motors  viewed  the  gas  in  the  combustion  chamber  when  the  charge 
is  burned  up  as  the  isentropic  gas  of  uniform  stagnation  and 
assumed  that  it  changes  according  to  the  isentropic  process 
(abbreviated  as  the  isentropic  model) .  However,  in  reality, 
the  gas  parameters  when  the  charge  is  burned  up  are:  the 
changes  of  the  velocity,  pressure,  temperature  and  density 
along  the  length  of  the  combustion  chamber;  because  of  the 
influence  of  the  added  flow  when  the  charge  burns,  the  entropy 
of  the  gas  on  each  section  of  the  combustion  chamber  in  the 
descent  segment  are  not  equal.  Therefore,  based  on  the  calcula¬ 
tion  of  the  isentropic  model,  how  large  are  the  errors  which 


can  finally  be  carried  along?  This  paper  uses  the  character¬ 
istic  method  to  consider  the  changes  of  the  gas  parameters 
along  the  length  of  the  combustion  chamber  to  calculate  the 
descent  segment  of  the  internal  ballistics  and  makes  compar¬ 
isons  with  the  calculations  of  the  isentropic  model. 

II.  Control  Equations 

Hypotheses : 

1.  We  only  consider  that  the  unsteady  flow  of  the  gas  in 
the  combustion  chamber  and  the  nozzle  flow  is  quasi-stationary. 
Because  the  nozzle  volume  is  much  smaller  than  the  volume  of 
the  combustion  chamber,  this  hypothesis  is  appropriate;  on  the 
other  hand,  in  order  to  reduce  the  computer  time,  we  consider 
that  the  unsteady  flow  of  the  nozzle's  changing  section  greatly 
reduces  the  length  of  the  computation  steps  and  causes  the 
computer  time  to  increase  several  times  or  several  tens  of  times. 

2.  The  combustion  gas  is  completely  gas  and  the  flow  is 
one-dimensional. 

3.  At  the  same  time  the  charge  on  each  section  of  the  com¬ 
bustion  chamber  is  burning  up,  this  hypothesis  makes  the 
computation  of  the  descent  segment  convenient.  If  we  have 
already  completely  computed  the  internal  ballistics  of  the 
motor's  operating  section  and  we  know  the  distribution  of  the 
gas  parameters  along  the  length  of  the  combustion  chamber  when 
the  charge  is  burned  up,  we  can  then  cancel  this  hypothesis. 

4.  The  combustion  chamber  is  a  uniform  gas  passage. 

5.  We  do  not  consider  that  the  combustion  of  the  adiabatic 
layer  is  equal  in  the  descent  segment. 

Under  these  hypotheses,  the  control  equations  are: 


Continuity  equation 


_*4L+._|P  +  p ±L  -o 
9t  9*  9* 


Momentum  equation 


Energy  equation 


Gas  state  equation 


Velocity  of  sound 


In  the  formulas,  u,  T,  P  ,  p,  a  and  y  are  the  velocity  of  flow, 
temperature,  density,  pressure,  velocity  of  sound  and  specific 
heat  ratio  of  the  gas.  This  is  a  set  of  hyperbolic  partial 
differential  equations  which  can  be  solved  by  the  characteristic 
method.  The  characteristic  equations  and  compatibility  equa¬ 
tions  are: 

Characteristic  equations 

(-£-). -*-7  <W) 

.  ,  (2) 
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Key:  (1)  Trajectory;  (2)  Mach  line. 
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Compatibility  equations 


4>.-«>4>.-0  (»**)  (1) 
dp*±psd*4~0  (fi4HI|)(2) 

Key:  (1)  Along  the  trajectory;  (2)  Along  the  Mach  line. 

In  the  formulas ,  corner  note  "0"  or  "t"  indicate  along  the 
trajectory  or  along  the  right  and  left  lines. 

This  paper  uses  the  contragradient  method  limited  difference 
algorithm  to  carry  out  computations  and  the  computations  of  the 
initial  values  are  determined  by  the  computation  of  the  added 
flow  the  instant  the  charge  burns  out. 

III.  Computation  Examples 

Fig.  1  shows  the  computed  motor.  Fig.  2  shows  the  calcula¬ 
tion  of  the  distribution  of  the  beginning  instantaneous  entropic 
increase  along  the  length  of  the  combustion  chamber.  This  is 
created  by  the  added  flow. 


L  L 


Fig.  1  Computed  motor* 
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Fig.  2  Distribution  of  entropic  increase  along  the  length  of 

the  combustion  chamber  when  t=0. 

Keys  (1)  Calories/kg* degrees;  (2)  Meters. 

Figs.  3a  and  3b  show  that  the  velocity  of  flow  of  gas 
changes  with  time.  Fig.  3a  shows  the  situation  of  the  descent 
segment  within  the  beginning  3  milliseconds  and  Fig.  3b  shows 
the  situation  within  4  milliseconds  after  60  milliseconds  in 
the  descent  segment.  We  can  see  from  the  figures  that  aside 
from  the  velocity  of  flow  of  the  front  end  being  zero  and  the 
rear  end  velocity  of  flow  monotonously  reduced,  the  velocities 
of  flow  on  the  other  sections  which  change  with  time  all  assume 
vibration  characteristics.  Because  the  times  the  disturbance 
waves  in  the  nozzle  entrance  area  are  propagated  to  each  sec¬ 
tion  are  different,  the  oscillation  characteristics  of  the 
velocity  of  flow  on  each  section  are  also  different.  We  can  see 
from  Fig.  3a  that  when  the  disturbance  is  propagated  to  a 
certain  section,  the  velocity  of  flow  of  that  certain  section 
has  maximum  value  and  after  the  disturbance  wave  has  passed,  its 
velocity  of  flow  then  decreases.  By  comparing  Figs.  3a  and  3b, 
we  can  see  that  the  oscillation  period  of  the  velocity  of  flow 
enlarges  with  the  passage  of  time  and  the  amplitude  of  the 
oscillation  decreases  with  the  passage  of  time.  We  can  see 
from  Fig.  3a  that  the  oscillation  period  is  1.8  milliseconds  and 
we  can  see  from  Fig.  3b  that  the  oscillation  period  is  2.15  mil¬ 
liseconds. 


Fig.  4  shows  that  the  distribution  of  the  pressure  and 
total  pressure  along  the  axis  of  the  combustion  chamber  change 
with  time.  Fig.  5  shows  that  the  average  pressure  of  the  com¬ 
bustion  chamber  changes  with  time.  We  can  see  from  Figs.  4 
and  6  that  the  general  tendency  of  these  parameters  which 
change  with  time  is  to  become  smaller  and  smaller.  However, 
they  are  distributed  along  the  axis  and  sometimes  the  front 
end  is  higher  than  the  rear  end  and  sometimes  the  rear  end  is 
higher  than  the  front  end  and  assumes  oscillation  character¬ 
istics.  As  shown  in  Fig.  5,  the  ratio  of  rear  end  pressure  p^ 
total  pressure  pgL  and  temperature  TL ,  and  front  end  pressure 
P1#  total  pressure  p^  and  temperature  changes  with  time,  its 
oscillation  period  and  velocity  of  flow  are  the  same  and  it  is 
the  time  the  disturbance  wave  is  propaged  back  and  forth  once 
in  the  combustion  chamber.  The  oscillation  period  within  the 
initial  several  milliseconds  of  the  descent  segment  is  1.8 
milliseconds. 


_ L,  ySL  and  _L  change  with  time. 


Key:  (1)  Milliseconds. 
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Fig.  6  The  combustion  chamber  pressure  changes  with  time. 

Key:  (1)  Calculated  according  to  the  idrnytopy; 

(2)  Calculated  according  to  the  characteristic 
method;  (3)  Milliseconds. 


Fig.  6  shows  the  average  combustion  chamber  pressure  calcul 
ated  according  to  the  characteristic  method  and  the  combustion 
chamber  pressure  which  is  calculated  according  to  the  isen- 
tropic  model  and  changes  with  time.  The  initial  value  of  the 


descent  segment  calculated  according  to  the  characteristic 
method  is  calculated  by  the  equilibrium  pressure  formula.  We 
can  see  from  the  figure  that  the  combustic  *  chamber  pressure 
calculated  according  to  the  isen tropic  model  is  slightly  lower 
than  the  value  calculated  by  the  characteristic  method.  Fig.  7 
shows  the  motor's  vacuum  thrust  taking  into  consideration  the 
unsteady  flow  in  the  combustion  chamber  and  the  motor's  vacuum 
thrust  which  is  calculated  according  to  the  isentropic  model 
and  changes  with  time.  Because  the  influence  of  the  unsteady 
process  in  the  combustion  chamber  on  the  thrust  is  very  small , 
the  thrust  curve  is  similar  to  the  combustion  chamber's  pressure 
curve  and  the  thrust  calculated  according  to  the  isentropy  is 
slightly  lower  than  the  value  calculated  by  the  characteristic 
method.  Calculations  show  that  the  differences  of  the  above 
curves  are  mainly  due  to  the  differences  of  the  initial  pressure 
(the  pressure  of  the  initial  point  of  the  descent  segment)  value. 
The  pressure  value  calculated  according  to  the  equilibrium 
pressure  formula  is  lower  than  the  average  combustion  chamber 
pressure  obtained  by  calculating  the  added  flow  in  the  uniform 
passage.  If  the  initial  combustion  chamber  pressure  calculated 
by  the  isentropic  model  is  equal  to  the  average  value  of  the 
initial  pressure  calculated  by  the  characteristic  method,  then 
the  calculation  results  of  the  two  methods  are  very  close. 
Therefore,  if  we  do  not  study  the  distribution  of  the  gas  para¬ 
meters  along  the  combustion  chamber  axis  but  only  study  the 
design  of  the  after-working  impulse  of  the  motor  or  combustion 
chamber  pressure  in  the  descent  segment  and  the  motor's  thrust 
which  changes  with  time,  we  can  calculate  according  to  the 
isentropic  model.  However,  the  initial  pressure  must  use  the 
average  combustion  chamber  pressure  obtained  from  the  added 
flow  in  the  uniform  passage. 
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Fig.  7  The  thrust  changes  with  time. 

Key:  (1)  Kg;  (2)  Calculated  according  to  the  isentropy; 

(3)  Calculated  according  to  the  characteristic  method; 

(4)  Milliseconds. 

IV.  Conclusions 

1.  Because  the  combustion  chamber  has  the  propagation  of 
disturbance  waves  originating  from  the  area  of  the  nozzle 
entrance,  in  the  combustion  chamber,  aside  from  the  front  end 

end,  the  gas  velocity  of  flow  on  each  section  which 
changes  with  time  assumes  oscillation  characteristics. 

2.  The  ratio  of  the  combustion  chamber's  rear  end  pressure, 
total  pressure  and  temperature,  and  the  front  end  pressure, 
total  pressure  and  temperature  which  changes  with  time  also 
assumes  oscillation  characteristics,  its  oscillation  period  is 
equal  to  the  time  required  for  the  disturbance  waves  to  be 
propagated  from  the  rear  end  to  the  front  end  and  then  from  the 
front  end  to  the  rear  end.  Because  the  propagation  speed  of 

the  disturbance  waves  in  the  descent  process  gradually  decreased, 
the  oscillation  period  enlarged  with  the  passage  of  time. 


3.  The  ranges  of  the  pressure  and  temperature  changes 
along  the  length  of  the  combustion  chamber  were  extremely 
small.  Moreover,  the  gas  velocity  of  flow  in  the  combustion 
chamber  was  very  low  and  its  influence  on  the  pressure  and 
temperature  was  very  small.  Therefore,  the  pressure  and  temp¬ 
erature  on  each  section  of  the  combustion  chamber  which  change 
with  time  do  not  resemble  the  oscillation  of  the  velocity  of 
flow  but  have  monotonous  descent.  However,  sometimes  the 
descent  is  relatively  fast  and  sometimes  it  is  relatively  slow. 

4.  If  we  do  not  study  the  gas  flow  parameters  in  the  com¬ 
bustion  chamber  which  are  distributed  along  the  axis  and  change 
with  time  but  only  consider  the  calculation  of  the  motor's 
after-working  impulse,  we  can  use  the  isentropic  model.  How¬ 
ever,  the  initial  pressure  must  use  the  average  pressure 
obtained  by  calculating  the  added  flow  when  the  charge  burns 
up. 

5.  Because  the  pressure  difference  of  the  combustion 
chamber's  front  end  and  rear  end  is  not  large,  the  gas  velocity 
of  flow  of  the  combustion  chamber's  rear  end  is  not  large,  the 
influence  of  the  unsteady  flow  in  the  combustion  chamber  on  the 
motor's  thrust  is  not  large  and  when  calculating  the  after¬ 
working  impulse  of  the  motor,  we  can  not  calculate  the  unsteady 
flow  in  the  combustion  chamber. 
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